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Numerical evaluation of the path integral for QCD on a discrete space-time lattice 
has been used to calculate ground state matrix elements specifying moments of quark 
density and spin distributions. This talk will explain how these matrix elements have been 
Q I calculated in full QCD using dynamical quarks, show how physical extrapolation to the 
chiral limit including the physics of the pion cloud resolves previous apparent conflicts with 
experiment, and describe the computational resources required for a definitive comparison 



O ■ with experiment. 
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1. INTRODUCTION 



The structure of hadrons differs profoundly from that of other familiar many body 
systems. Unlike electrons in atoms or nucleons in nuclei, quarks are confined in the 

^ ■ nucleon. Whereas photons can be subsumed into the Coulomb potential in atomic physics 
and pions can be subsumed into the nucleon-nucleon potential in nuclear physics, gluons 

k> ! are essential degrees of freedom in light hadrons that carry half the momentum and 
H ' have important nonperturbative topological excitations. Since nonperturbative QCD is 
presently intractable analytically, the goal of this work is to use lattice field theory to solve 
QCD with controlled errors to provide a quantitative understanding of the rich quark and 
gluon structure of the nucleon and to obtain insight into how QCD actually works to 
produce this structure. 

Decades of high energy scattering experiments have now provided detailed experimen- 
tal knowledge of the light cone distributions of quarks and gluons in the nucleon, so this 
work addresses the use of lattice QCD to understand the observed distribution of the 
quark density and helicity. Using the operator product expansion, it is possible to cal- 
culate moments of quark distributions, and I will discuss here the first calculations in 
full QCD|l|Q]. A major puzzle in the field has been the fact that quenched calculations 
of these moments, which ignore quark-antiquark excitations of the Dirac sea, disagree 
with experiment at the 20-50% level. I will show that contrary to some conjectures, at 
the quark masses accessible in practical calculations, including quark loops does not alter 
the results significantly. Rather, I will argue that the physical origin of the discrepancy 
with experiment has been incorrect extrapolation to the physical quark mass, and will 
show how extrapolation incorporating the leading non-analytic behavior required by chiral 
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Figure 1. Connected (upper row) and disconnected (lower row) diagrams contributing to 
hadron matrix elements. The left column shows typical contributions of quarks and the 
right column shows contributions of antiquarks 

symmetry produces consistent results for the moments of quark distributions. In addi- 
tion, we have also compared full QCD results with configurations that have been cooled 
to remove all the gluon contributions except for those of instantons and shown that the 
qualitative behavior of the moments is reproduced by the instanton content of the gluon 
configurations. 

2. MOMENTS OF QUARK DISTRIBUTIONS IN THE PROTON 

By the operator product expansion, moments of the following linear combinations of 
quark and antiquark distributions in the proton 

' da;x'^(g(a;) + (-l)"+^g(x)) 






(1) 



dxx"(Ag(x) + (-l)"Ag(x)) 



dxx"(5g(x) + (-l)"+^5g(x)) 



are related to the following matrix elements of twist-2 operators 

{PS\tlj-f^^HD^^---iD^"^\PS) = 2{x''-^)gP^^'---P^-^ 
(P5|?/^7^^i75zL>^2 ■ • • iD^'"^\PS) = 2(x""^)Ag MS^^'^P''^ ■ ■ ■ P'^"^ 
(P5|^at"^^^l75zP)^^ ■ • ■ iD^'-^\PS) = 2{x''-^)iq MS^^'P^^'^P"^ ■ ■ ■ P^-\ 



(2) 



Here, q = q\ + qi, Ag = q-[ — qi, 6q = qT + q±, x denotes the momentum fraction carried by 
the quark, and { } and [ ] denote symmetrization and antisymmetrization respectively. We 
note that odd moments (a;")q are obtained from deep inelastic electron or muon scattering 
structure functions Pi or P2, even moments of {x"')^q are determined from gi, and these 
moments are proportional to the quantities f „+i and a„ defined in ref . . In addition, g2 
also determines the quantity dn 



{PS\'il:-f^''-f5iD^^''^ ■ ■■iD''-^\PS) = -dnS^'^P^^''^ . . . p^- 

n 



(3) 




Figure 2. Overlap between smeared source and proton ground state as a function of the 
source RMS radius. The overlap for zero smearing is approximately 10"'^ 

Even moments {x^)q are obtained from deep inelastic neutrino scattering, and in addition, 
a variety of other processes have contributed to what is now a detailed empirical knowledge 
of the quark and antiquark distributions in the nucleon. We will subsequently compare 
our results with moments calculated from the CTEQ, MRS, and GRV global fits to the 
world supply of data. 



3. CALCULATION OF MATRIX ELEMENTS 

Proton matrix elements of the operators in Eq.^ are calculated by evaluating the con- 
nected and disconnected diagrams shown in Fig. |l]. Note that both the connected and 
disconnected diagrams each receive contributions from quarks and antiquarks. Depend- 
ing on the moment, by Eq.^, the sum of the diagrams yields either the sum or difference 
of the moments of the quark and antiquark distributions. In the past, there has been 
some confusion on this point: because the connected diagrams are called the valence 
quark distribution and experimentalists define the valence quark distribution as the dif- 
ference between the quark and antiquark distributions, the connected contributions have 
erroneously been compared only with the difference between the empirical quark and 
antiquark distributions. Because it is technically much more difficult to evaluate the 
disconnected diagrams, our present calculations only include connected diagrams. For- 
tunately, the disconnected diagrams are flavor independent, so they cancel out of the 
difference between up and down quark distributions. Hence, in Table |^, we compare lat- 
tice calculation of the difference between connected diagrams for up and down quarks with 
the corresponding difference of moments of experimental data for the sum or difference 
of quark and antiquark distributions. 

On the lattice, connected diagrams are evaluated by calculating a three point function 
in which a source creates a state with the quantum numbers of the proton, the operator 
acts on this state, and a sink flnally annihilates the state. Because evolution in imaginary 
time filters out the ground state, when the operator is sufficiently far from both the source 
and sink, it acts in the ground state and produces the desired ground state matrix element. 
As the time at which it acts approaches either the source or sink, it sees excited state 
contaminants, yielding a central plateau corresponding to the physical matrix element 
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Figure 3. Plateaus in measurements of operators in a zero momentum ground state as a 
function of the euclidean time separation from the source. 

and exponential contaminants at the edges. Obviously, it is beneficial to optimize the 
overlap of the source with the ground state to maximize the plateau region and minimize 
the effect of the excited state contaminants at the edges. 

In this work, connected diagrams were calculated using sequential propagators gener- 
ated by the upper two components of the nucleon source J" = u'^%{C'^^)i3^pid^ e"'^'^ . The 
overlap with the physical proton ground state was optimized using Wuppertal smearing 
[^] to maximize the overlap -P(O) = |(J|0)p. Figure Q shows that varying the smearing 
reduced P by over 4 orders of magnitude, yielding an overlap with the physical ground 
state of approximately 50%. Dirichlet boundary conditions were used for quarks in the 
t-direction. 

The resulting plateaus for four operators that could be measured in a proton with 
zero three-momentum are shown in Fig. ^. Here one observes both a statistically well 
determined central plateau region and the effectiveness with which the excited state con- 
taminants have been reduced by the optimized source. 



4. OPERATORS AND PERTURBATIVE RENORMALIZATION 

The continuum operators defined above are approximated on a discrete cartesian lattice 
using representations of the hypercubic group that have been chosen to eliminate operator 
mixing as much as possible and to minimize the number of non-zero components of the 
nucleon momentum. The operators we have used are shown in Table 1, where we have 
indicated whether the spatial momentum components are non-zero and whether mixing 
occurs. Note, no* indicates a case in which mixing could exist in general but vanishes 
perturbatively for Wilson or overlap fermions and no** indicates perturbative mixing 
with lower dimension operators for Wilson fermions but no mixing for overlap fermions. 
Because the statistical errors are much larger for sources projected to non-zero momentum, 
the moments corresponding to operators requiring non-zero momentum are presently not 
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Figure 4. Comparison of linear chiral extrapolations of full and quenched calculations of 
(xq) and (xAq) showing agreement within statistical errors 



well determined. 

To convert from lattice regularization at the scale of the inverse lattice spacing 1/a 
to the continuum MS scheme at momentum scale Q, we use the one-loop perturbation 
theory result 

The anomalous dimensions 7,^ and the finite constants Bij we have calculated and used 
in this work are tabulated in Table 2 0. 



Table 1 Lattice Operators 
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Table 2 Perturbative renormalization 

constants 
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Figure 5. Comparison of linear chiral extrapolations of full and cooled calculations of (xq) 
showing the extent to which instantons reproduce the full result. 
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Figure 6. Chiral extrapolation of {x)u~d using Eq. ^ 

5. RESULTS 

The moments listed in Table 1 were calculated [|l| on 16'^ x 32 lattices for Wilson 
fermions in full QCD at /5 = 5.6 using 200 SESAM configurations at each of 4 k's and 
at /5 = 5.5 using 100 SCRI configurations at 3 k's. They were also calculated with two 
sets of 100 full QCD configurations cooled with 50 cooling steps and in quenched QCD at 
(3 = 6.0 using 200 configurations at each of 3 k's. Typical linear chiral extrapolations for 
operators calculated with nucleon momentum equal to zero are shown in Figure ^ for full 
and quenched calculations of {x)q and {x)Aq, showing agreement within statistical errors. 
To avoid finite volume errors at the lightest quark mass, the SESAM [H results were 
extrapolated using the three heaviest quark masses. Table 3 shows a major result of our 
work, that there is complete agreement within statistics between full and quenched results. 
Statistics with the SCRI configurations 0] are not yet adequate to present extrapolations 
in the coupling constant. 

Typical chiral extrapolations for cooled configurations are compared with the corre- 
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sponding uncooled full QCD calculations in Figure |^. This qualitative agreement be- 
tween cooled and uncooled results occurs at light quark mass for all the twist-2 matrix 
elements we calculated and demonstrates the degree to which the instanton content of 
the configurations and their associated zero modes dominate light hadron structure Pj. 

A longstanding puzzle in this field has been the fact that when quenched lattice results 
are linearly extrapolated in rrig, results disagree at the 20-50% level. Our results show 
that inclusion of quark-antiquark excitation from the Dirac Sea does not resolve this 
discrepancy. Salient examples from Table |^ are {x)u-d ~ 0.25 — 0.27 compared with the 
experimental result 0.15 and qa = (l)Au-Ad ~ 1-0 — 1.1 compared with the experimental 
result 1.26. 

There is strong evidence that the physical origin of these discrepancies is the inadequate 
treatment of the pion cloud in the nucleon that has been necessary because of limited 
computational resources. By necessity, present calculations are restricted to quark masses 
that are so heavy that the pion mass is above 600 MeV and the pion cloud surrounding 
the nucleon is strongly suppressed. Physical quantities like the nucleon magnetic moment 
and axial charge clearly depend strongly on the pion current, and should therefore be 
very sensitive to the absence of the full pion cloud. Furthermore, because of the rapid, 
nonlinear variation from the chiral logs arising from Goldstone boson loops, it is clear 
that a linear extrapolation is completely inadequate to describe the correct chiral physics. 

In a recent work[10|, also discussed by A. Thomas at this conference, we have shown 
that chiral extrapolation incorporating the leading non-analytic behavior from chiral per- 
turbation theory can systematically resolve the discrepancies in the moments {x"')u-d 
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using the formula: 



(3(7^2 + 1)^2 



2 I 1^^2 ^2 



Inf- "^^ 



+ &„m2 (4) 



(47r/^)2 Vrn2+^2, 

The coefficient of the leading non-analytic behavior m1.ln{m\) is determined from chiral 
perturbation theory. The parameter /i specifies the scale above which pion loops no longer 
produce rapid variation. It corresponds to the upper limit of the momentum integration if 
one applies a sharp cutoff in the loop integral and physically corresponds to the inverse size 
of the quark core of the nucleon that serves as the source for the pion field. As shown in ref 
r^ , the value /i ~ 550 MeV, which is consistent with the value required to extrapolate 



the nucleon magnetic moment and with chiral nucleon models, extrapolates the world's 
supply of lattice data to the experimental values of {x)u-d, {x'^)u-d, and {x^)u-d- In fig ^ 
we show the extrapolation of our lattice data for {x)u-d- From this figure, it is clear that 
5% measurements down to m^ = 0.05 GeV^, would provide data for a definitive lattice 
calculation. This calculation will require 8 Terafiops-years and thus can be carried out on 
the next generation of 10-Teraflops computers. 
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